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We report on a work in progress, whose goal is a systematic field theoretical deriva- 
tion of the quantum transport equations for baryon production in the electroweak 
plasma at a first order phase transition in the limit of slowly varying background 
^ ' fields (thick wall limit). We start with the Schwingcr-Dyson equations for the two 

Cf^ , point Green function written in the closed time contour (CTC) formalism. The 

quantum Boltzmann equations for the density matrix arise when the SD-equations 
^^ ■ are expanded to the first order in the gradients in the on-shell limit. In this paper 

\y^ ' we consider only scalar fields, but the formalism can easily be extended to fermions. 

o' 

0^ 



1 Introduction 



,S^ ' The matter-antimatter asymmetry to be explained is most accurately pinned 

Mh, down by nucleosynthesis considerations.j-AUowing for a fair uncertainty in the 

O ' present observational situation, one hasOJ 

(D : 

^.- ,75^ — =.2-7x10-". (1) 
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An exciting realisation is that rjB could have been generated at the electroweak 
phase transition. What is required is a strong transition and an adequate 
C^ ' source of CP-violation.j^ While the minimal standard model (MSM) is known 

to fail these conditions □, theyp.may well be met in its popular minimal super- 
symmetric extension (MSSM)H. There are three distinct parts of the problem: 
(1) equilibrium calculations, which provide information about the strength of 
the phase transition El; (2) phase trcmsition dynamics, whose study provides 
the phase boundary shape and speedu; and finally (3) baryon production com- 
putation. The goal of our program is to bring problem (3) to the level of 
understanding of (1). 

At the moment there is no systematic quantum field theoretical derivation 
of the CP-violating sources and the transport equations that would consistently 
describe all aspects of the problem in one formalism. The difficulty lies in the 
inherent out-of-equlibrium nature of the problem. The theory we study here, 



is defined by the following scalar lagrangian 

L = {d^<fy\d^cf) - 0tm2(a;)0 - A„t , (2) 

where >Cint contains the interaction terms. The field can have several compo- 
nents so that rri} is in general a complex matrix varying in space and time, pos- 
sibly giving rise to out-of-equilibrium conditions and nontrivial CP-properties, 
which are the necessary ingredients for baryogenesis. 

2 The Kadanoff-Baym equations 

We start our analysis with the exact SchssdiCiger-Dyson equations in the Keldysh 
closed time contour (CTC) formalism laoQ. The basic quantity is the out-of- 
equilibrium 2-point Green function 



Gf{x,V) = -ilTc\K{Mfi{v^ (3) 



where Tc defines time ordering along the contour C which starts at some io 
(often taken to be at — oo) goes to -|-oo, and then back to ^q. From now on 
we suppress the indices a, /3, which denote different scalar fields ("flavours"), 
and work with matrices in this flavour space. The two point function Gc{x, y) 
obeys the contour Schwinger-Dyson equation 

Gc{x,y) = Gl{x,y)+ f dx' j dx" GI{x,x')^c{x\x")Gc{x" ,y) , (4) 
Jc Jc 

where Ec is the self-energy and Gq is the free particle (tree-level) propagator. 
In order to solve Gc{x, y) from (0), some additional information about Ec niust 
be provided. In general Ec is specified through the usual BBGKY hierarchy 
of equations for n-point functions. In the weak coupling limit Ec can be 
computed perturbatively, if the interaction lagrangian Cj in Eqn. (0) is given; 
for the moment however, we can keep Ec completely general. 

After some algebrja Eqn. (0) can be recast as the well-known Kadanoff- 
Baym (KB) equationsu: 

(Go 1 - E'-''^) (g) C'"^ = S 
(Go ^ - E'-) (g) G<^> = E<'> (g) G"'. (5) 

where (E> denotes the convolution integral: [/ ® g]{x, y) = J dx' f{x, x')g{x' , y) 
and G'"'° are the retarded and advanced Green functions and 

G>{x,y) = -i{cp{x)cj,^{y)), G<{x,y) = -i{4>\y)cp{x)) (6) 

are the quantum Wigner functions. The Kadanoff-Baym equations (ph are still 
exact, but formidable to solve. They considerably simplify when expanded to 
the first order in gradients. 



3 Gradient expansion 

Gradient expansion is valid in the limit of slowly varying background fields. 
This means that the Compton wave length £ ^ k^^ oi excitations should be 
smaller than the characteristic length scale of the background variation. More 
precisely, propagators and self-energies should vary slowly with respect to the 
macroscopic (average) coordinate X — {x + y)/2, i.e. 

\G\ » \dxdkG\ , m » \dxdk^ , (7) 

where the precise meaning of momentum k is defined below. Applying this 
gradient expansion to Eqns. (ph we get the new KB-equations 



e 
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{Go^-E'-^'^llG"'''^} = 1 (8) 

e-^^{G-i - J:^}{G<^>} = e-**{l]<^>}{G"^} , (9) 

where the Green functions G = G{k;X) and self-energies S = T,{k;X) are 
now defined in the mixed representation, defined by the Fourier transform 
with respect to the relative coordinates r ^ x — y (Wigner transform) : 

G{k; X)^ I d^re'''''G{X + r/2, X - r/2) . (10) 

The O-operator of Eqn. (0) is defined as the following generalization of the 
Poisson brackets: 

0{f}{9} - \ [dxf ■ dkg - d^f ■ dxg] ■ (H) 

To the first order in gradients, e"** in Eqn. (^) simplifies to 1 — i<>. For the 
electroweak phase boundary the typical macroscopic scale is of the order of the 
boundary thickness L ^ (10 — 20)/tBu, whereas for a typical plasma excitation 
dk ^ k^^ ^ T^^, so that dxdk ^ l/LT ^ 1, i.e. criteria (Q) are satisfied. 

To the lowest order in gradients equations (jg^ clearly define the fa- 
miliar propagators of an excitation in a translationally invariant background, 
and their poles define the well known quasiparticle dispersion relations. Thus 
Eqn. (pi) in particular contains information about the spectrum of excitations 
of the system, and it is contained in the spectral function defined by 

A^'-{G'--Gn='^{G>-G<). (12) 



4 Propagator equation 

We begin by considering the propagator equation for the retarded and advanced 
Green functions. For the one-field case Eqns. (||) imply 

cosO{Go-i-E'''"}{G''^'^} = 1 (f3) 

sinO{Go"i-S''''^}{G''^'^} =0. (f4) 

Now observe that Eqn. (^^ can be obtained from Eqn. (^3|) by an applica- 
tion of the differential operator tanO, and hence it gives no new information. 
Therefore, to the first order in gradients, we have simply 

G"'" ^{G^^ -T."''")-^. (f5) 

From this, one can easily read off the spectral function A (|lj) and Gr — 
(G'- + G'»)/2. 

In the case of several flavours Eqn. ( [l5|) does not solve the propagator 
equation (^ to the first order in gradients. For simplicity consider S^ = 
(S'' -1- S'')/2 -^ 0. Then the explicit solution has the following form 

(^r.a^^t ^^^ [^] K (IQ) 

det[X]\ + k,[[H„\]_,\]_±ieuj{det[\]+XtiiX]) ' ^ ' 

where U is the unitary matrix in the flavour space that diagonalizes the mass 
matrix to^, H^ — —iUdzU^ is a "gauge" field, and A = U{Gq^ — ^r)U^ is 
the diagonal matrix with the entries \± = w'^ — uj'i^ , o;^ = A;^ + m\ , and m\ 
axe the diagonal entries of m\, — Unn?U^ for the two field case. Note that 
an inhomogeneous background alters the pole structure of the propagator. 
Indeed, Eqn. (^6|) can be diagonalized, and one can show that each of the 
homogeneous particle poles (jJ± is split into three poles, each with a weight 
1/3 of the original excitation, which lie on the circle in the complex frequency 
plane, with shifts given by {i5w,(5we^''^/^,(5we~^*'^/^}, where 5uj oc {Hz)i2 is a 
nonanalytic shift. One can show, however that these nonanalytic shifts affect 
the on-shell Boltamann equation at the second order in gradients and hence 
can be neglectedll3. Finally we note that it is not hard to generalize Eqn. (16) 
and include the self-energy T,r. 

5 On-shell approximation for the Quantum Bolzmann Equation 

The hermitean part of Eqn. (0) defines the quantum Boltzmann equation 
(QBE), and anti- hermitean part a nondynamical constraint equation (CE). 



The QBE is the dynamical master equation whose solutions specify how the ef- 
fective phase space of the system is populated, and the CE is the non-dynamical 
constraint which singles out the physical solutions. In the on-shell limit the 
QBE (upper signs) and OSE (||) simphfy to 

-tO{G^' - Sfi}{G<}^ + {Go 1 - Sh}{G<}^ 

- Ki^"}iG<}± - {^<}{G>} J + {S<}{GhL (17) 



where {a}{&} — ab^ ha. In the following we only consider the QBE. 

Eqns. ( |l7|) specify the dynamics of the quantum Wigner function G^. In 
the on-shell limit the singular structure of G< can in general be factored out 
as follows 

iG< = A"- + nMs = 7rsign(w)^g'(5(w-Wi), for T -> , (18) 

where n stands for the (regular) matrix of "occupation numbers." 

Upon integrating the QBE from Eqn. dT^) over positive frequencies, we 
obtain the following on-shell quantum Boltzmann equation for particles: 



1 1 i 

■^{dj:Cjd,dsfd}- ■^{dsoJd,df:fd} + ^, 



dtf + 7;{dk'^d, dsfd} - T^jS^cDd, 9,-/4 + -[^d, fd] + Rot fd = Coll[fd]. (19) 



This is our main result. Analogous equation for antiparticles is obtained by 
integrating over negative frequencies. Here ivd is the (diagonal) matrix of the 
quasiparticle frequencies. The derivative terms, including the anticommuta- 
tors, form the hermitian many-field generalization of the usual Boltzmann flow 
term, and the commutator term is the expected mixing term familiar from the 
usual Liouville equation. This term is the source of the nontrivial quantum 
coherence effects, and together with some additional rotational terms included 
in Rot f it leads to potential new sources for baryogenesis. The collision term 
Coll[fd] in ( p9| ) depends on the detailed form of the interactions. In addition 
to the usual diagonal collision integrals it contains new terms which tend to 
damp the off-diagonal elements of fd to zero. 

Due to the nonlocal character of the SD equation, fd is related to C^'^ in 
a nontrivial manner in that there is a derivative correction: 

/d = ^ {'^d,5d} - ^ [h^d, dsQd] -f — [ds(^d: d^gd] , (20) 



where gd = lAgW is defined in Eqn. (18). This definition is the unique one 



that renders the QBE in the simple form of Eqn. (pi 



The basic structure of Eqn. ( ]l9| ) is basis independent; when rotated into 
another basis by an unitary matrix U, fd and Cjd become / =Li'^ fdti and Cj = 
U^'cudU, respectively, where w in general is not diagonal. Further, the detailed 
form of TZot term changes. Physically the TZot term arises as a consequence of 
the space-time dependent rotation required to bring an excitation to the local 
WKB basis; discussion of its detailed form is beyond the scope of this talk. 

In the adiabatic limit, where the wall width L is much larger than the in- 
teraction rates, L ^ l/F, collisions efficiently damp away the off-diagonal ele- 
ments (unless interactions are flavour blind), resulting in quantum decoherence. 
In this case the quantum terms vanish, and one obtains the diagonal semiclas- 
sical Boltzmann equations for the^quasiparticle distribution functions, which 
were studied earlier in refs.EJ andEI, with only the classical force Fd — —d^LOk 
and "spontaneous" terms to fuel baryogenesis. However, we conclude by not- 
ing that even when TL ^ 1 is satisfied, it may be important to include the 
non-adiabatic effects, when we are interested in phenomena that vanish in the 
adiabatic limit, as may be the case with CP violating effects; the result is then 
suppressed by (TL)^^. 

We thank Dietrich Bodeker, Felipe Freire, Guy Moore, Emil Mottola, 
Michael Schmidt, and Misha Shaposnikov for many discussions and comments. 
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